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Chiral transition and deconfinement in QCD 
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The study of QCD with two light dynamical fermions is of fundamental importance to understand the mech- 
anism of color confinement. We present results of a numerical investigation on the order of the chiral phase 
transition with Nf = 2 by use of a novel strategy in finite size scaling analysis. We compare the critical behaviour 
of the specific heat, of the chiral susceptibility and of the equation of state with the possible critical behaviours. 
A second order transition in the 0(4) and 0(2) universality classes are excluded by our data and substantial 
evidence emerges for a first order transition. Like in most of previous works we have used the standard staggered 
action with L± = 4: possible scaling violations and the need for further studies are discussed. 
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1. Introduction 

The study of the deconfinement transition 
in presence of two light degenerate dynamical 
quarks (Nf = 2) is of special interest: the sys- 
tem is a good approximation to Nature and it is 
very interesting also from a theoretical point of 
view. A sketch of the phase diagram for Nf = 2 
is shown in Fig. ^ m is the quark mass and [i is 
the baryon chemical potential. 

In the /i = plane, quarks decouple in the 
limit m — * oo and the system tends to the 
quenched limit, where the deconfining transition 
is an order-disorder first order phase transition; 
Z 3 is an associated symmetry and the Polyakov 
line (L) is an order parameter. The inclusion of 
dynamical quarks explicitely breaks the Z3 sym- 
metry and (L) is not a good order parameter, 
even if it works as such for quarks masses down 
to m ~ 2.5 - 3 GeV. 

At to ~ a chiral phase transition takes place 
at a critical temperature T c ~ 170 MeV, from a 
low temperature phase in which chiral symmetry 
is spontaneously broken to a high temperature 
phase in which it is restored: the corresponding 
order parameter in this case is the chiral conden- 
sate (ijjil>). Also the ?7a(1) symmetry, which is 
explicitely broken by the axial anomaly, is ex- 
pected to be effectively restored at some tempera- 



* Speaker at the Conference 



ture T4 >T C . Empirically also the Polyakov line 
has a rapid increase at the same transition tem- 
perature as for chiral symmetry, indicating de- 
confinement: it is not clear which relation exists 
between the chiral transition and the deconfining 
transition. The transition line depicted in Fig. ^ 
is defined by the maxima of a number of suscep- 
tibilities (Cy, Xrn, ■ ■ ■), which indicate a rapid 
variation of the corresponding parameters across 
the line: the positions of these maxima all coin- 
cide within errors. 

At to ~ it is possible to perform a renor- 
malization group analysis plus e-expansion tech- 
niques, assuming that the relevant degrees of free- 
dom for the chiral transition are scalar and pseu- 
doscalar fields |1I2I3| . The result is that if the 
Ua(1) symmetry is effectively restored, i.e. if the 
77' mass vanishes at T c , then there is no IR sta- 
ble fixed point and the phase transition is first 
order; if this is not the case an IR fixed point 
exists, which can produce a second order phase 
transition in the 0(4) universality class. 

If the first case is realized, the transition is 
first order also at to ^ and most likely up 
to to = 00. In the second case a phase transi- 
tion is present only at m = 0, while a continuous 
crossover takes place as to ^ 0: that means that 
one can move continuously from confined to de- 
confined and that no true order parameter exists. 
However this would be in contradiction with the 
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experimental knowledge we have about confine- 
ment: indeed the upper limit on the number of 
free quarks per proton is R = n q /n p < 10~ 27 , 
while R ~ 10~ 12 is expected from the Standard 
Cosmological Model. A reduction factor 10~ 15 is 
difficult to explain in natural ways unless R = 0, 
which means that confinement is related to some 
symmetry of the QCD vacuum and therefore the 
deconfining transition is associated to a change 
of symmetry of the vacuum, i.e. it is an order- 
disorder phase transition rather than a continu- 
ous crossover. 

The determination of the order of the chi- 
ral transition for Nf = 2 is therefore a funda- 
mental issue, which can be settled by numeri- 
cal simulations of lattice QCD. The problem has 
been investigated by several groups with stag- 
gered |4I5I6I7I8I9I1()| and Wilson fermions. 
The strategy used has been either to look for signs 
of discontinuity at the transition, or to study the 
scaling with respect to m of different susceptibil- 
ities and of the magnetic equation of state. No 
clear discontinuities have been observed, but also 
no conclusive agreement of scaling with 0(4) crit- 
ical indexes. A general tendency exists however 
in the community to consider the chiral transition 
second order, and the line of Fig. ^ a crossover. 
We present the results of a big numerical effort 
using large lattices, made in order to clarify the 
issue. A full account of our results can be found 
in Like most of the other works we use non 
improved Kogut-Susskind action, lattices 4 x 




Figure 1. Schematic phase diagram of Nf = 2 
QCD. 



with L s = 12, 1 6, 20, 24, 32, and the non exact R 
hybrid algorithm for Nf = 2 Scaling viola- 
tions are expected and a more careful study with 
L t — 6, with an improved action and algorithm 
is planned in order to control them. 

2. Results 

The theoretical tool to investigate the order of a 
phase transition is finite size scaling. The extrap- 
olation from finite size L s to the thermodynam- 
ical limit is governed by critical indexes, which 
identify the order and the universality class of 
the transition. 

Approaching the transition, for a higher or- 
der or weak first order transition, the correlation 
length of the order parameter £ goes large com- 
pared to the lattice spacing a, so that the depen- 
dence of physical quantities on a/£ can be ne- 
glected. More precisely, if C/kT is the effective 
action (density of free energy) 
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the dependence on a/£ disappears as T c is ap- 
proached, since £ diverges as 



£ ~ T ^oo (2) 

where r = 1 — Sr. Therefore around the chi- 
ral transition the system has two fundamental 
lengths: the correlation length £ and the inverse 
quark mass l/m q . The variable L s /^ can be 

traded with Tl}J v and the scaling law follows 

^cL;U(rLl^,am q L^) . (3) 

As r — > irrelevant terms can be neglected and 
the correlators of the order parameter describe 
the thermodynamics. The most important quan- 
tity is however the specific heat, which shows the 
correct critical behaviour independently of the 
identification of the order parameter. 

The following scaling laws hold for the the spe- 
cific heat and for the susceptibility of the order 
parameter 



Cv-C'o^ L«/»<t> c {tL]/ v 



, am q L y s h j , 



, am q L v s h 



(4) 
(5) 
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Co steins from an additive renormalization. It is 
also possible to write them in the alternative form 

C v - Co ~ Lf^ (r(am q )- 1 ^ vy ''\am q L^ (6) 

X ^LZ' v ^ x [r{am q )- 1 /^\am q L^) . (7) 

The phase transition is characterized by the 
critical indexes: the values relevant to our anal- 
ysis are listed in Tabled 0(4) is the symmetry 
expected if the transition is second order, How- 
ever lattice discretization can break it down to 
0(2) for Kogut-Susskind fermions at non zero 
lattice spacing. 

Finite size scaling analysis is made difficult by 
the presence of two independent scales. In the 
previous literature the assumption has been usu- 
ally made that the volume is large enough so to 
neglect the dependence on L s : since at fixed am q , 
[3 the susceptibilities must be analytic in the ther- 
modynamical limit, as L s goes large the depen- 
dence on am q L v s h must cancel the dependence on 
L 8 in front of the scaling functions in Eqs. (J2J and 
@. It follows that 

C v - C ~ (am q )- a /^f c (rlam,)" 1 /^)) (8) 

x ~ (om,)-^)/ x (rH,)- 1 /^)) . (9) 

The peak values of (Cy — Co) and of x should 
then scale as 

(Cy - C ) max oc (am,)-"/^' 

Xmax oc (am,)^) (10) 

as am q — > 0. The pseudocritical couplings, i.e. 
the positions of the maxima, should instead scale 
as 

T(am q )- 1/( - uyh) = const. (11) 
Table 1 
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One can also consider to keep tLJ fixed while 
taking aL s 3> l/m w . This assumption should 
work better if L s is still comparable to the corre- 
lation length, which may be the case close enough 
to the critical point. In this case the scaling laws 
are 

C v - Co ~ K)-°/W/ c (tLII») (12) 

x ~ (am q )-^^f x (tL 1 /") . (13) 

Eqs. (|1L)|) stay unchanged, the positions of the 
maxima scale as 

tL\I v = const (14) 

while the widths of the peaks in this case are vol- 
ume dependent. 

We have instead followed a novel strategy 
which does not rely on any assumption: in or- 
der to reduce the problem to one scale we have 
kept fixed one of the scaling variables and we 
have studied the dependence on the other. As 
one can see in Tabic ^ the index y h is the same 
for 0(4) and 0(2) symmetry. We have therefore 
made a number of simulations at different values 
of L s and am q keeping am q L y s h fixed and assum- 
ing yh = 2.49 which corresponds to 0(4) or 0(2). 
In this way as L s is increased, am q — > 0, so that 
the infinite volume limit corresponds to the chiral 
transition at am q — 0. 

From Eqs. (J2J and (JHJ it follows that the max- 
ima at constant am q L v s h scale as 

(Cy-CoW oc L a J v 

Xmax OC L?J V . (15) 

as L s — > oo and their positions scale as 

tL\I v = const (16) 

If 0(4) or 0(2) is the correct symmetry, the val- 
ues of a/v and 7/1/ should be consistent with the 
corresponding values listed in Tabled 

We have run two such sets of Monte Carlo sim- 
ulations, called in the following Runl and Run2, 
with am q L y s h = 74.7 and am q L v s h = 149.4 respec- 
tively. The spatial lattice sizes L s used for each of 
the two sets are L s = 12, 16, 20, 32, the standard 
hybrid R algorithm has been used to update 
configurations. 



4 



0.3 

B ,« 0.25 
U° 0.2 
0.15 
0.1 



1 


i 1 i 


II 


1 




• 0(4| scaling 








■ ■ 0(2| scaling 




A - 




Run 1 












z ~ 


























T 




ii 


1,1,1 



> 0.15 

8 



- 1 1 


1 1 1 


II 


i 




• 0(4| scaling 








■ ■ 0(2| scaling 




A - 




Run 2 












l - 


























1 









12 14 16 



20 22 24 26 28 30 32 

L 



20 22 24 26 28 30 32 

L 



Figure 2. Specific heat peak value for Runl (left) and for Run2 (right), divided by the appropriate powers 
of L s to give a constant. Both the 0(4) and 0(2) critical behaviors are displayed. 
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Figure 3. The same as figure for the chiral susceptibility Xn 



In Figs. [5] and [3] we show the peak values of 
the specific heat and of the chiral susceptibility, 
divided by the power of L s appropriate for 0(4) 
or 0(2), as a function of L s (see Ea. 1151 : scal- 
ing is clearly violated, 0(4) and 0(2) universality 
classes are excluded by our data. 

An alternative way to investigate the order of a 
phase transition is to study the scaling of pseudo- 
critical couplings: one can try the two alternative 
scaling laws of Eq. (fHEjl . i.e. t c — k T L 1 J v , or of 
Eq. (TU, i-e- t c = ^(om,) 1 /^) . One has to be 
careful in defining the reduced temperature r in 
presence of dynamical fermions, since in this case 
the physical temperature T = \/{L t a{j3, m q )) is a 



function of both (3 and am q , so that the reduced 
temperature r can be expanded as a power series 
in (/3 — Pq) and in am q , where /3q is the chiral 
critical coupling. Only the linear term in (3 was 
considered in the previous literature. We have 
found that the following terms are sufficient to fit 
the data 

r cx (0Q— P) + k m am q + 

+ k m 2(am q ) 2 + k m pam q (f3 ~ 0) ■ (17) 

As a result of our analysis we have found that it 
is not possible, within the present mass range, to 
discriminate among the various possible critical 
behaviours by looking at pseudocritical couplings 
only. We would like to remark that the inclusion 
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of other terms in Eq. (|17|l . besides the one linear 
in P solely considered in previous literature, is 
crucial to obtain any scaling at all. 

We have also tried to make a scaling analysis 
in the same way as done in previous literature, 
i.e. supposing that the lattice size is much larger 
that all other relevant physical lengths. No uni- 
versality class is chosen a priori in this case and 
one can test all the possible critical behaviours. 
We have found that assuming that aL s 3> 1/TOtt 
but still L s — £, i.e. Eqs. JT2J), lO works bet- 
ter: this is reasonable around the critical point, 
where £ goes large. We have added to the data 
from Runl and Run2, those from two other simu- 
lations performed at L s — 16, am q — 0.01335 and 
L s = 24, am q = 0.04444. In Fig0|we show the 
scaling obtained for the specific heat peak: 0(4) 
is again clearly excluded, while a good agreement 
is found with a weak first order critical behaviour. 
An analogous behaviour is observed for the chiral 
susceptibility. 

As a further test of scaling we check the equa- 
tion of state 

{tpip) ~ m 1 ' 5 f{T{am q y l K vy ^) (18) 

No scaling whatsoever is observed, neither 0(4)- 
0(2) nor first order, if the raw measured data 
are introduced in Eq. © . After a proper subtrac- 
tion ^21 we have instead studied the following 
scaling law 

(^) - (H>)o = (am^FiTiam,)- 1 /^) (19) 



Results are shown in Fig. again the first or- 
der behaviour describes well the data while the 
second order is excluded. 

3. Conclusions 

We have argued that the study of the order of 
the chiral phase transition for Nf — 2 is of fun- 
damental importance to understand confinement. 
By adopting a novel finite size scaling strategy 
which reduces the analysis to a one scale problem, 
we have been able to show clear incompatibility 
of our data with 0(4) (0(2)) second order critical 
behaviour. After introducing the correct defini- 
tion of the reduced temperature, we have shown 
that the analysis of pseudocritical coupling alone 
cannot discern between the possible critical be- 
haviours. 

We have also repeated the same finite size scal- 
ing analysis performed in previous literature, i.e. 
assuming that the thermodynamical limit has 
been reached, and again we have found disagree- 
ment with 0(4) (0(2)) second order critical be- 
haviour and consistency with a weak first order 
critical behaviour both in the scaling of suscepti- 
bilities and in the equation of state. This would 
be in agreement with confinement being an ab- 
solute property of matter related to some sym- 
metry and with the deconfining transition being 
order-disorder, even if we have still not found any 
clear evidence for discontinuities in physical ob- 
servables. 
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Figure 5. Scaling of the equation of state for 0(4) (left) and first order (right). 



We would like to remark that our simulations, 
like in most of previous works, have been per- 
formed with the standard staggered action at 
L t = 4 and with the non exact Hybrid R algo- 
rithm: scaling corrections could be important in 
this case and the issue is surely still open: we plan 
in the future to investigate it more deeply by mak- 
ing simulations with improved actions and algo- 
rithms and with am q L v s h fixed according to first 
order. 
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